Valence-Bond Dynamical Mean-Field Theory of Doped Mott Insulators 
with Nodal/ Antinodal Differentiation 
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We introduce a valence-bond dynamical mean-field theory of doped Mott insulators. It is based 
on a minimal cluster of two orbitals, each associated with a different region of momentum space and 
hybridized to a self-consistent bath. The low-doping regime is characterized by singlet formation 
and the suppression of quasiparticles in the antinodal regions, leading to the formation of Fermi 
arcs. This is described in terms of an orbital-selective transition in reciprocal space. The calculated 
tunneling and photoemission spectra are consistent with the phenomenology of the normal state of 
cuprates. We derive a low-energy description of these effects using a generalization of the slave-boson 
method. 

PACS numbers: 71.27.+a,71.30.+h,74.72.-h 
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The doping of a Mott insulator is a fundamental 
problem of condensed matter physics, relevant to the 
physics of cuprate superconductors [lj. In the simplest 
Brinkman-Rice Q description, the doped metallic state 
is a Fermi liquid in which quasiparticles are formed with 
a heavy mass m*/m ~ 1/8 and a reduced weight Z ~ 5 
(6 is the doping level). This physical picture can indeed 
be rationalized using the modern theoretical framework 
of dynamical mean- field theory (DMFT) It is ap- 

plicable when spatial correlations are weak, which is fa- 
vored by high dimensionality and strong competing (e.g. 
orbital) fluctuations. In cuprates however, which are two- 
dimensional materials with low orbital degeneracy, it was 
pointed long ago by Anderson in a seminal paper [l[ 
that the antiferromagnetic superexchange (J) plays a 
key role, leading to strong short-range correlations as- 
sociated with singlet formation (valence bonds) between 
nearest-neighbor lattice sites. Slave-boson mean-field 
theories [M, la, H, Bl , as well as projected variational wave- 
functions [9j, llO( > provide simple theoretical frameworks 
to incorporate this effect, modifying the Brinkman-Rice 
picture at small doping S < J/t and leading, e.g. to 
a finite effective mass m*/m ~ l/(J/t + S), consistent 
with observations in cuprates. However, these theo- 
ries fail to describe a key phenomenon in underdoped 
cuprates, namely the strong differentiation in momen- 
tum space observed e.g. by photoemission spectroscopy 
(ARPES) 11,]: Coherent quasiparticle excitations are 
suppressed in the antinodal regions of the Brillouin zone 
(BZ) and a pseudogap appears in the normal state. In 
order to take this phenomenon into account while incor- 
porating short-range correlations, cluster extensions of 
the DMFT framework have been investigated by several 
groups 0, [H, EH • Most studies have considered clusters 
of at least four sites (plaquette) and numerical efforts 
have been devoted to increase the cluster size in order 



to improve momentum- resolution and get closer to the 
two-dimensional lattice 
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In this article, we follow a different route, looking for a 
description based on the minimal cluster able to success- 
fully describe momentum-space differentiation together 
with Mott physics. We find that a two-site cluster is 
sufficient to achieve this goal on a qualitative level, and 
to a wide extent on a quantitative level when compared 
to larger cluster calculations. This allows us to con- 
struct a valencc-bond dynamical mean-field theory (VB- 
DMFT) of nodal/antinodal differentiation, in which this 
phenomenon is linked to the distinct properties of the 
orbitals associated with different regions of momentum- 
space. 

The main motivation to choose the smallest possi- 
ble cluster is to advance our qualitative understand- 
ing. Since the theory is based on a two-site Anderson 
model, results can be interpreted in terms of valence- 
bond singlet formation and linked to the competition 
between singlet-formation and individual Kondo screen- 
ing 
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However, the self-consistency of 
the bath does bring in novel aspects to this competition 
in comparison to the non self-consistent two-impurity 
model. The present VB-DMFT can be viewed as an ex- 
tension of the static mean-field theories (e.g. slave-boson 
based) of singlet formation. In contrast to those theories 
which have a limited number of static variational param- 
eters, it involves a dimer coupled to a self-consistent bath 
through energy-dependent hybridization functions. This 
allows for a description of the physics over a wide range 
of energy scales. At low energy however, the new slave 
boson approximation introduced in (l9j reproduces with 
a remarkable accuracy several aspects of the full solution. 

We study the Hubbard model on a square lattice, with 
hopping between nearest-neighbor (t) and next-nearest- 
neighbor sites (£'). In the following, we use U/t = 10 



2 



and t 1 /t = —0.3, which are values commonly used for 
modeling hole-doped cuprates in a single-band frame- 
work. All energies (and temperatures) are expressed 
in units of D = 4t = 1, and the doping is denoted 
by 8. We use a two-site effective Anderson impurity 
problem, involving the on-site interaction U and two hy- 
bridization functions: a local one An(w) = A22(w) and 
an inter-site one Ai2(u;), which are self-consistently de- 
termined by relating the two-impurity problem to the 
original lattice one. We have investigated several such 
embeddings, both of the dynamical cluster approxima- 
tion (DCA) and cellular-DMFT (CDMFT) type @, EH 
with similar results. Here, we focus on a somewhat gen- 
eralized form of the DCA embedding, which preserves 
the symmetries of the square lattice, in which the Bril- 
louin zone is decomposed into two patches of equal sur- 
face: a central square (denoted P+) centered at momen- 
tum (0,0) and the complementary region (P_) extend- 
ing to the edge of the BZ and containing in particular 
the (7r,7r) momentum. From the lattice Green's func- 
tion, two coarse-grained Green's functions in momen- 
tum space are constructed: G±{u>) — X^keP ± ^O 4 , w ) 
(with momentum summations normalized to unity within 
each patch). Following the DCA construction, the in- 
ner (resp. outer) patch self-energy is associated with the 
even- (resp. odd) parity self-energy of the two-impurity 
effective problem, i.e to the even (resp. odd) orbital com- 
binations (c\ ±cl)/y/2. Indeed, the states close to (0,0) 
have more bonding character while those close to (tt, tt) 
have more antibonding character. The self-consistency 
condition reads: Gk{u) = Y^k.£P K — ^ k (<*>)] • 
In this expression, the index K — ± refers both to the in- 
ner/outer patch index and to the even/odd orbital combi- 
nations. We solve the self-consistent two-impurity prob- 
lem using both continuous-time quantum Monte Carlo 
(CTQMC) 20] which sums the perturbation theory in 
A a b(iu! n ) on the Matsubara axis, and an approximate 
method geared at low-energy properties: the rotation- 
all y in variant slave-boson formalism (MSB) presented 
in [19| . The MSB method introduces slave-boson ampli- 
tudes 4>r n , a density matrix connecting the eigenstates 
|r) of the isolated dimer to the quasiparticle Fock states 
\n), determined by minimizing (numerically) an energy 
functional. 

In Fig. [TJ we display the real part of the even- and 
odd-orbital self-energy at zero frequency, as determined 
by both methods, as a function of 8. We find a rather re- 
markable agreement between the CTQMC solution and 
the low-energy MSB. The two orbitals behave in a simi- 
lar way at high doping 8 > 25%. Below this doping level, 
we observe an onset of orbital differentiation, which is a 
manifestation of momentum differentiation in the lattice 
model. This differentiation increases as 8 is reduced, un- 
til a transition is reached at 8 ~ 16% (in CTQMC). At 
this characteristic doping, fj, — £1(0) reaches the band 
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FIG. 1: (Color online) Left: real part of E± (0) as a function of 
doping level, computed with RISB (lines) and CTQMC (sym- 
bols), fi — Si(0) (diamonds) reaches the odd-orbital band 
edge (dotted line), which becomes empty at low energy be- 
low S ~ 16%. Right: statistical weights of the various dimer 
cluster eigenstates. S is the intra-dimer singlet, 1+ the (spin- 
degenerate) state with one electron in the even orbital, E the 
empty state and T the intra-dimer triplet. f3 = 100. 

edge corresponding to the odd orbital, and the latter be- 
comes empty at low energy and remains so for all lower 
dopings. G(k, u) no longer has poles at lo = in the 
outer patch, and low-energy quasiparticles exist only in- 
side the inner patch. Hence, at low doping, momentum- 
space differentiation becomes strong and manifests itself 
as an orbital-selective transition in VB-DMFT. 

In order to gain further qualitative insight, we also 
plot in Fig. [T] (right part) the statistical weight of several 
cluster eigenstates given within slave bosons by the 
amplitude pr = J2n \4>rn\ 2 - We compare it to a simi- 
lar estimate [2l[ from CTQMC. The agreement between 
CTQMC and MSB is again very good, and even quanti- 
tative for the two states with highest weights. At large 
doping, the empty state and the two spin-degenerate 
states with one electron in the even orbital dominate, as 
expected. As doping decreases, these states lose weight 
and the intra-dimer singlet prevails, reflecting the strong 
tendency to valence-bond formation. The states with 
immediately lower weights are the one-electron states 
and the valence-bond breaking triplet excitation which 
dominates over the empty state. Therefore, the orbital 
(momentum) differentiation at low doping is governed by 
intra-dimer singlet formation, reminiscent of the singlet 
regime of the two-impurity Anderson model. 

The gaping of the odd orbital (outer patch) is actu- 
ally a crude description of the pseudogap phenomenon. 
To illustrate this, we compute the tunneling conduc- 
tance dl/dV as a function of voltage V. This calcula- 
tion is made possible by the high quality, low-noise, of 
the CTQMC results on the Matsubara axis, allowing for 
reliable analytical continuations to the real axis at low 
and moderate energy, using simple Pade approximants. 
The conductance is displayed on Fig. [2] together with 
the gap A in the odd Green's function, obtained from 
A = Si(A) + e min — (j,, with e m i n the lower edge of the 
band dispersion in the outer patch. Note the overall 
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site cluster (plaquette), using the standard DC A embed- 
ding (with the BZ divided into 4 patches centered around 
K= (0,0),(0,7r),(7r,0),(7r,7r)). The results of Fig. [3] re- 



FIG. 2: (Color online) Left panel: tunneling spectra for dif- 
ferent doping levels. All curves are obtained using Eq. (6) of 
Ref. [22|] with the same proportionality factor. Right panel: 
gap A to the band edge of the unoccupied odd-orbital states. 
P = 100. 



particle- hole asymmetry of dl/dV and the peak at posi- 
tive voltage. This peak shifts towards higher energy with 
decreasing doping, as does the gap A, and can indeed be 
traced back to the edge of the unoccupied odd-orbital 
spectral function. These observations compare favorably 
to tunneling experiments in the normal state of under- 
doped cuprates [jjj ]. 

We now address two related issues: how to reconstruct 
information in momentum space using our two-orbital 
description, and how to gauge the reliability of a descrip- 
tion based on only two momentum-space components, 
as compared to calculations with larger cluster sizes and 
better momentum-space resolution. The approximation 
of lattice quantities from the cluster ones is a central issue 
in cluster methods. Periodization is crucial in CDMFT 
to restore translational invariance. In DCA, translational 
invariance is not broken, but there is still a large freedom 
when interpolating the self-energy in the BZ from the 
cluster self-energies. The most local quantity is expected 
to give the more reliable interpolation. We investigate 
two possibilities among those that have been discussed 
in the literature [23|: i) interpolating the self-energy (E- 
interpolation) as E(k, uj) — £ + (cj)a+(k) + E_(w)a_(k), 
with a±(k) = i{l ± i;{cos(k x ) + cos(fc y )]}; ii) inter- 
polating the cumulant (M-interpolation) as M (k, us) = 
M + (w)a + (k) + M_(w)a_(k). The cumulant is related 
to the self-energy by E = to + /i — M — 1 . It is the dual 
quantity of the self-energy in an expansion around the 
atomic limit and a natural measure of how much the 
hybridization to the self-consistent environment changes 
the impurity Green's function as compared to an iso- 
lated dimer. Close to the Mott insulator, it is more local 
than the self-energy and a better quantity to interpo- 
late [2^]. Comparing the interpolations on small clus- 
ters with direct calculations on larger clusters, having 
better k-resolution, provides a systematic test of clus- 
ter schemes and interpolations. As a first step, we com- 
pare in Fig. [3] the results of VB-DMFT, with S-orM- 
interpolation, to the cluster components Sjf(u) of a four- 




FIG. 3: (Color online) Comparison of the plaquette (dots) 
and dimer self-energies at K = (0, 0), (71", it), (tt, 0). Upper 
and lower-left panels: ImEif(i0 + ) as a function of doping. 
Lower-right panel: E„-o(ik->n) at fixed 5 = 0.08. The dimer re- 
sults at K = (71-, 0) are obtained by M-interpolation (squares) 
and E-interpolation (diamonds). At K = (0, 0), (n, n) both 
interpolations coincide. /3=100. 

veal two main points: i) the M-interpolation of the two- 
orbital results is clearly superior to the E-interpolation 
for reconstructing plaquette cluster quantities and ii) 
when M-interpolated, the two-orbital description does 
quite a remarkable job at capturing the full frequency- 
dependence of the various cluster components E^-(w) of 
the plaquette results. Note that the plaquette cluster- 
momentum K = (7r, 0) is not present as an individual 
orbital in the two-site description: it is entirely recon- 
structed by interpolation, and as such is the most direct 
test of the reconstructed momentum-dependence. A dis- 
tinctive feature of the results depicted in Fig. [3] is that 
the scattering rate near momentum (it, 0), ImE w o(*0 + ), 
displays a maximum around a doping level S ~ 8%, as 
previously noted in the plaquette study of [24[. Note 
however, that this maximum does not induce a maximum 
of the scattering rate computed at the Fermi surface. 

VB-DMFT provides a simple description of momentum 
differentiation as observed in ARPES experiments. This 
is illustrated by the intensity maps of the spectral func- 
tion A(k, 0) (obtained with M-interpolation) displayed in 
Fig. [U At very high doping S > 25% (not shown), clus- 
ter corrections to DMFT are negligible and the spectral 
intensity is uniform along the Fermi surface. In contrast, 
at lower 5, momentum differentiation sets in, revealing 
apparent "Fermi arcs" at finite temperature with higher 
spectral inten sity in the nodal direction in comparison 
to antinodes [H, [H, EE H3- The last panel in Fig. g] 



shows that the contrast of the spectral intensity along 
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FIG. 4: (Color online) Intensity maps of the spectral func- 
tion ^4(k, 0) for different doping levels. Lower-right panel: 
normalized intensity A(<j>, Q)/A(0, 0) along the Fermi surface 
(4> — is the node, <j> = ±45 the antinode). The nodal inten- 
sity A(0,0) is 0.045 for 5=6%, 1.66 for 5=10% and 4.61 for 
6=14%. (3 = 200. 



as a minimal cluster-based description of momentum- 
space differentiation in doped Mott insulators. Because 
of its simplicity, this theory can be investigated with 
moderate numerical effort and progress in qualitative un- 
derstanding can be achieved with low-energy methods 
such as rotationally invariant slave bosons. The cal- 
culated STM and ARPES spectra are consistent with 
the phenomenology of the normal state of cuprates. 
The low-doping regime is dominated by singlet forma- 
tion. Mott physics is responsible for the suppression 
of coherent quasiparticles at the antinodes, in quali- 
tative agreement with other approaches starting from 
the weak/intermediate coupling viewpoint (3^ |. Within 
VB-DMFT, this suppression is described as an orbital- 
selective transition in momentum-space. 
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the Fermi surface has a maximum around S w 10%, simi- 
larly to ARPES experiments (cf. Fig. 3B of \2J§). At low 
doping, singlet formation induces a large real part in 
(cf. Fig. [1]) and a large imaginary part of the self-energy 
in the (tt, 0) and (tt, it) regions, which are responsible for 
this strong momentum-space differentiation. At interme- 
diate doping (10% < 5 < 20%), this differentiation is reli- 
ably addressed using VB-DMFT. At low doping (5 < 8%) 
the M-interpolated self-energy develops singularities on 
lines in momentum space, leading to lines of zeroes of the 
Green's function and to the breakup of the Fermi surface 
into pockets [H S, El EH. A better momentum 
resolution (larger clusters) is necessary to obtain reliable 
results in this regime. 

VB-DMFT and the (non-self-consistent) two-impurity 
Anderson model share common features. In both cases, 
at low-6, the singlet state dominates, and the real part 
of the odd-orbital self-energy is large. These effects are 
due to the term transferring singlet pairs from the even 
orbital to the odd orbital, as can be checked by explic- 
itly removing it from the dimer Hamiltonian. Interest- 
ingly, strong fluctuations in the singlet pairing channel 
and momentum-space differentiation appear to be re- 
lated effects. The key difference between VB-DMFT and 
the two-impurity model with fixed bath is that the self- 
consistency leads to the opening of a gap in the odd or- 
bital. This gap reduces the scattering rate of the even 
orbital, leading to a maximum in ImS_|_(«0 + ) (and also 
in the reconstructed ImS w o(*0 + )), which is absent in the 
non-self-consistent two-impurity model. 

To summarize, we have proposed in this article a 
valence-bond dynamical mean-field theory (VB-DMFT) 
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